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Definition 1. Let A be a square, nonsingular matrix. ‘Then 
the inverse matriz A~! of A is a unique matrix for which, 


AAMT =A- A 
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Example 1. An inverse matrix may be found using the for- 
mula, 
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Example 2. Matrix equations of the form 
Ax=b 
can be solved with the help of the inverse matrix A7! as 
x=A'b 


where A is an n X n matrix, x a vector of size n whose com- 
ponents are variables, and b a vector of size n containing con- 
stants. 


Business mathematics, Linear algebra, 22" November 2005 —3- From 5” November 2005 , as of 10°” December, 2005 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Theorem 1. Let A be the coefficient matrix and A; a matrix 
formed from A by replacing the column of coefficients of x; with 
the column vector of constants. Cramer’s rule solves a system 
of linear equations through the use of determinants as follows. 
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Definition 2. Let a system of n functions not necessarily 
linear be 


ViL= fi Misesagen) 


Ui pO igasts en) 
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Then a Jacobian determinant comprises all the first-order par- 
tial derivatives of the system arranged in ordered sequence, that 


1S 


Oy 
Bae E OY)... ia = ea 
OX),.. OX OYn 


Ox} 
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Theorem 2. Let a system of nm equations be 


ST Disanceea) 
(ea Cera 7 


If | J| =0, then y; are functionally dependent. 
On the other hand if | J| #0, then y; are functionally indepen- 


dent. 
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Definition 3. A determinant | H | composed of all the second- 
order partial derivatives, with the direct partials on the prin- 
cipal diagonal and the cross partials off the same, is called a 
Hessian. In other words, let a multivariable function be 


z= f(x,y) 
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Then the Hessian of z is 


Z Z 
| H | _ | “£2 LY 
eyt “yy 
where Zr, = 2,7. Moreover, the first principal minor is 
y y Pp P 


dy, | = 2e5 


and the second principal minor is 


| He | = _ i = Zeezyy — (Zay)” 
ry Zyy 
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Theorem 3. Let a multivariable function be 
2= Joy) 
and let the first-order conditions 
are met. Then a sufficient condition for z to be at optimum is 
2 
ZenZyy > (Zay) 
together with 
Zea, Zyy <0 
in case of a maximum and 
Zea, Zyy > OV 
in case of a minimum. 
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Definition 4. From Definition 3, 


if | H,| > 0 and | Ha| > 0 the Hessian | H | is said to be positive 
definite, and the second-order conditions for the minimum are 
met. 


If | H,| <0 and | H2| > 0 it is said to be negative definite, and 
the second-order conditions for the maximum are met. 
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Algorithm 1 Procedure to test for the optimality of multivariable functions 
of two variables. 


c= f(x, y) 
find z, and z, 
if z, =O and z, = 0 then 
find Z,,, 2, and Z,, 
find H, and H, 
if | H,| > 0 and | H,| > 0 then 
|H| is positive definite 
elseif | H,| <0 and | H, > 0| then 
| H| is negative definite 
endif 
endif 
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Definition 5. Let y = f(21,...,%,) be function of n vari- 
ables. Then the n"-order Hessian for this function is 


Yl o¢':) Yin 
|H|= cee ; 
YUnl °*° YUnn 


Then the first principal minor | Hy | is simply x1;, and the 2” 
principal minor is 


Yiloccs Yili 


Yil cc° Via 
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Theorem 4. Let y = f(x1,...,2n) be function of n variables. 
Let the Hessian of y be represented by | H |. 


Then if all the principal minors of | H | are positive, then | H | is 
positive definite and the second-order conditions for a relative 
minimum are met. 

If the sign of the principal minors alternates between negagive 
and positive, then | H | is negative definite and the second-order 
conditions for a relative maximum are met. 
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Example 3. For 
y = f (x1, £2, £3) 
the third-order Hessian is 
Y11 4 Y12 = Y13 


|H|=|]ye1 yoo 23 
Y31  Y32 —-Y33 


where 
Oy Ay 


Ww ==, ¥12 = ———, and so on 
Ors @ 029024 
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The first-, second- and third-order Hessian’s are respectively 


Y1l Y12 


ie Hy | = 
|| =y11, | He | os 


and 
Yl Yi2 = -Y/i13 
|H3| =] yo1 Yyo2 23 
Y31 Y32 Y33 
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If |H,| > 0, |H2| > 0 and | H3| > 0, then H is positive 
definite and the second-order condition for minimum is fulfilled. 


If | H,| < 0, | H2| > 0 and | H3| < 0, then | H | is negative def- 
inite and the second-order condition for maximum is satisfied. 
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Definition 6. <A discriminant is a determinant of a quadratic 
form. Let the quadratic form be 


z= ar’ + bxy + cy” 


which is in matrix form 


z=[x y| 


Nic & 
oO NIMS 
8 


Then the discriminant is 


= 


Nico & 
D bvloa 
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The first principal minor of the discriminant is 


| Di |=a 
and the second principal minor 


| D2| = 


Nico & 

oO NI& 
o 
ND 
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Theorem 5. Let a quadratic form be 
z= ar’ + bxy + cy” 


and let the discriminant of z be | D|. 


If |D,| > 0 and | D2| > 0, then | D| is positive definite and 
z > 0 for all x,y £0. 


If |D,| < 0 and | Dg| > 0, then | D| is negative definite and 
z <0 forall z,y 40. 
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Theorem 6. _ Let 


fay) 


be a function subject to a constraint 


g(x,y) =k 


where k is a constant. ‘hen the optimisation of f can be done 
by first transforming f together with g into a new function 


F(z,y,r) = flz,y) +A(k— 9(@,y)) 


Business mathematics, Linear algebra, 22" November 2005 —21-— From 5” November 2005 , as of 10°” December, 2005 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


and then solve the following equations, 


F(a, Y; X) 
Fy (a, Y, x) 
Fy(a, y, >) 


| 
SSO 


to obtain the critical values x9, yo and Ag at which F’ and hence 
f are optimised. 
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Definition 7. In the constrained optimisation with Lagrange 
multipliers in Theorem 6 above, f is called an 


objective or origin function 


and F' the 


Lagrangian function 


Business mathematics, Linear algebra, 22" November 2005 —23- From 5” November 2005 , as of 10°” December, 2005 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 8. Let 
TGisntag 


be a function of n variables subject to constraints 


CEs Da) 
Let 


FE Vissers Dyn) = 7 eat eg tn A =H 6 insist ny) 
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Then the bordered Hessian | H | is defined as either 


fi, Fy 
7 F591 
|| = 


Or 
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This is simply the Hessian 
py. “eee die, 


F Gt. P= Finn 
bordered by the first derivatives of the constraint with zero 
on the principal diagonal. The order of a bordered principal 
minor being determined by the order of the principal minor 
being bordered, 7 7 
|| = | A | 
since in this case an n X n principal minor is being bordered. 
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Theorem 7. Let f(21,...,%) be a function of n variables 
subject to constraints g(11,...,2n). 


Let | H| be the bordered Hessian defined in Definition 8. 


Then if | Hy|,...,|Hn| < 0, then the bordered Hessian | H | 
is positive definite, and therefore is a sufficient condition for a 
minimum. 


If | H2| > 0 | H3| < 0, | H4| > 0, and so alternatingly on, then 
|H| is negative definite, which is a sufficient condition for a 
maximum. 
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Example 4. Let f(z, y) be a function to be optimised subject 
to a constraint g(x,y) = k, where k is a constant. Then the 
Lagrangian function becomes 


F(z,y,r) = flz,y) +A(k— 9(@,y)) 


The first-order conditions for optimisation are 


F, =F, =F, =0 
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The second-order conditions for optimisation can be expressed 
together as a bordered Hessian 


a Peg Foy Gx 
| HT | = Pyar Fyy Gy 
Gx Dy 0 
or 


_ 0 Gx Dy 
| 7 | = gx Fre Foy 


Jy Fyx Fyy 
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Note 1. ‘Theorem 6 gives the first-order conditions for opti- 
mising a function subject to some constraints. ‘Theorem 7 gives 
the second-order conditions for optimising a function subject to 
some constraints. 
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Definition 9. A Marshallian demand function gives an ex- 
pression of the amount of a good that a consumer will buy as a 
function of commodity prices and income available. It is derived 
by maximising the utility subjected to a budgetary constraint. 
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Example 5. Let a utility be 
U = 4142 
which is subject to a constraint 


Pigi + p2gq2 = 6 


where 6 is the amount of income available, that is to say, our 
budget. Then the Lagrangian function is 


U = qqo t+ A(b— pig — p2q2) 
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The first partial derivatives are then 
uy = gg — Ap, = 0 (1) 
ug = q — Ap2 = 0 ( 
uy = b— pig — p2g2 = 0 ( 


where uy, ug are respectively ug, and Ugo. 
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Simultaneously solving Equation’s 1, 2 and 3 leads us to 


P\ P2 


Hence go = q1p1/p2 and q, = qop2/p; and from Equation 3 we 
have, 
b= pig + pe = poe + DI 
P2 Pi 
which yield us, for gj and qo, the Marshallian demand functions 
which maximise satisfaction of the consumer subject to income 
and prices. 
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Next, we test the second-order conditions by firstly finding 
uj1 = 0, ua2 = 0, Wig = var = 1, 91 = pi and go = pe, which 
give us 


_ 0 Ll py 
\H#|=|1 O po 
pi p2 O 


which gives | H)| = 2p;p2 > 0 Hence | | is negative definite 
and thus u is maximised. 
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Definition 10. ‘The production process of producing one good 
usually requires the input of many other intermediate goods. 
Let x; be the total demand for product 7, and let b be the final 
demand for the product from the ultimate users. Then, 


C;,= ajyjxry~t...+ Ajntn + 0; 


for 2 = 1,...,n, where aj; is a technical coefficient which rep- 
resents the value of input 7 required to produce one monetary 
unit’s worth of product 7. 
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If we consider the total demand for every one of the products, 
then 
x= Ax+b 


where 
zi G11 o<'' Gin by 
In Gnl "** Gann On, 
It follows from this that 
x =(I—A)'b 
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The matrix A is known as the matrix of technical coefficients. 
It is also known as the input-output table, the rows being the 
inputs and the columns the outputs.. The matrix J/—A is known 
as the Leontief matriz. 
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Example 6. In a complete input-output table, labour and 
capital would also be included as inputs. ‘These give the value 
added by the firm. They are normally put as an extra row at the 
bottom of the matrix of technical coefficients A. The vertical 
summation of each column of the table is then equal to 1. 
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Definition 11. Let A be a square matrix. Then a scalar 
such that the equation 


Av = Av (4) 


holds for some vector v ~ O is called an eigenvalue + of A, 
and the vector v is called an ezgenvector of A corresponding 
to the eigenvalue A. The eigenvalue A is also known as the 
characteristic root, or the latent root, while the eigenvector is 
also known as the characteristic vector, or the latent vector. 
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Note 2. From Equation 4 it follows directly that 
(A —Al)v =0 (5) 


Then A — Al is called the characteristic matrix of A. Since v 
assumes a unique value and by assumption v + 0, it follows 
that A — AI must be singular, which means that its rows must 
be a multiple of one another. Now A — AI is zero if and only if 
the characteristic determinant | A — I | of A is zero. 
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In other words 
| A—AI| =0 (6) 
which is called the characteristic equation of A. With Equation 
6 there will be an infinite solution for v in Equation 5. In 
particular, if v is a solution, that is if it is an eigenvector, so 
is kv for any k # 0. We force a unique solution by using the 


normalisation 
Ee 
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Then the sign-definiteness of A can be determined from the 
characteristic roots A’s. 


Thus if all \’s are positive, then A is positive definite; and if 
negative, negative definite. 


Let at least one A be zero, which is neither positive nor nega- 
tive, if all the remaining A’s are nonnegative, then A is positive 
semidefinite; and if they are nonpositive, negative semidefinite. 


Lastly, if some of the \’s are positive while others are negative, 
then A is indefinite. 


Business mathematics, Linear algebra, 22" November 2005 —43-— From 5” November 2005 , as of 10°” December, 2005 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Note 3. We have seen in Note 2 how, having found 4;, where 
1=1,...,n, we find through normalisation the corresponding, 
unique v;. On the other hand if we have found first the v;,’s, 
their corresponding A;’s may be found by first forming a trans- 
formation matriz 
Pe ea Vp 

and then the corresponding eigenvalues or the characteristic 
roots are obtained from 


1 OO -:: 0 
T' AT =| - a 
0) oe 


Business mathematics, Linear algebra, 22" November 2005 —44— From 5” November 2005 , as of 10°” December, 2005 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 12. The vector equation, Equation 4, has as its 
solutions the zero vector v = 0 together with all the correspond- 
ing elgenvalue-eigenvector pairs. ‘The set of all the eigenvalues 
of A is called the spectrum of A. The spectral radius of A is 
then the largest of all the absolute values of the eigenvalues of 
A, that is to say, 


max |); 
t 


The set of all eigenvectors v;;, together with 0, forms a vector 
space called the ezgenspace of A corresponding to 4;. 
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